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Abstract 

In this paper we use the matrix string approach to begin a study of high energy 
scattering processes in M-theory. In particular we exhibit an instanton-type configuration 
in 1+1 super- Yang-Mills theory that can be interpreted as a non-perturbative description 
of a string interaction. This solution is used to describe high energy processes with non- 
zero longitudinal momentum exchange, in which an arbitrary number of eigenvalues get 
transferred between the two scattering states. We describe a direct correspondence between 
these semi-classical SYM configurations and the Gross-Mende saddle points. We also study 
in detail the pair production of D-particles via a one-loop calculation which in the 1+1D 
gauge theory language is described by the (perturbative) transition between states with 
different electric flux. Finally, we discuss a possible connection between these calculations 
in which D-particle production gives important corrections to the Gross-Mende process. 
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1. Introduction and Summary of Results 

High energy processes in string theory were first considered from the point of view of 
conventional string perturbation theory by Gross and Mende [U in the regime of fixed angle 
scattering and in the near forward regime by Amati et al in 0. The recent insights from 
M-theory, however, have provided a large number of new non-perturbative tools which 
can now be used to put these works into a new perspective, and extend the results into 
new directions. For instance, it was long believed that the string length l s marks the 
minimal distance that can be probed via scattering processes in string theory. This belief 
was based on the fact that fundamental strings tend to increase in size when boosted to 
high energies, and thus appear to be incapable of penetrating substringy distance scales. 
Since the discovery of D-particles as non-perturbative solitons of the IIA theory, however, 
we know that there exists small scale structure that, at least for weak string coupling, 
extends well below the string length This particular realization provided important 

motivation for the Matrix theory conjecture of 00 that all localized excitations of M- 
theory (including the fundamental strings) are representable as multi-D-particle bound 
states 0,0. 

In this paper we begin a study of high energy processes in type IIA string theory, 
by making use of this Matrix theory formalism. We focus on the four graviton scattering 
amplitude, and in particular we will present a detailed calculation of the pair production 
rate of D-particles via this process. Our aim is to probe in this way the transition re- 
gion between the conventional perturbative string regime and the strong coupling regime 
described by 11-dimensional M-theory. (See fig. 1.) 

From the ten dimensional perspective of IIA string theory, D-pair production is an 
inelastic scattering process, in which two strings exchange one unit of D-particle charge. It 
is inherently nonperturbative and thus inaccessible to conventional perturbative methods. 
It is also inaccessible in the traditional Matrix theory approach since the anti-D particles 
are boosted to infinite energy. From the eleven dimensional perspective, on the other 
hand, the D-pair creation process can simply be thought of as the elastic scattering of two 
particles in which one unit of Kaluza-Klein momentum in the 11 direction is exchanged. 
Via this interpretation, one can rather straightforwardly obtain a tree level estimate of 
the probability amplitude. This estimate should be reliable for large values for the S 1 
compactification radius Rn and for collision energies sufficiently below the 11-dimensional 



1 For reviews see [f -10] and references therein. 
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Planck energy. At high energies and/or small values for -Rn , on the other hand, we expect 
the physics of the scattering process to be quite different from (semi-) classical supergravity. 
In the following we will attempt to gain more insight into this regime via the Matrix string 
approach. 
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Fig. 1: The phase diagram of the S 1 compactification of M-theory, with horizontal 
axis the log of the length scale and vertical axis the log of the string coupling. The 
various perturbative and low energy limits are indicated. The shaded region marks 
the regime where D-pair production is expected to be the dominant high energy 
process. 



Matrix string theory arises from the original Matrix theory proposal [[7] via compact- 
ification on a circle, and starts from the action of 1+1-dimensional maximally supersym- 
metric Yang-Mills theory with gauge group U(N), 

S= [dr fdaTrl-^-F", - \- (D^X) 2 + ^[X,X? 

J Jo k 4 z w s n i\ 

+i$ P1> - -jT* [X it ^]\ , 
9s > 

defined on the circle < a < 1. The derivation of this action is reviewed in the Appendix; 
here and henceforth we work in string units, l 8 = 1. Via the identification of the eigenvalues 
of the matrices X with the transverse location of type IIA supersymmetric strings, this 
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SYM model can be reinterpreted as a non-perturbative formulation of light-cone gauge 
IIA string |T3| - |T5[1 . In this correspondence, the string coupling constant g s is inversely 
proportional to the Yang-Mills coupling g YM and the free string limit therefore arises in 
the strong coupling limit of the Yang-Mills model. This correspondence has been worked 
out in some detail in |15| . More generally, however, all regimes of the S 1 compactification 
of M-theory, as indicated in fig 1, should according to the Matrix string conjecture of 
|p], p^| , p^JT5| 1 via the above identifications be described by particular regimes of (the large 
iV limit) of the 1+1D supersymmetric gauge theory. 

In particular, it is expected that in the weak coupling, moderate energy limit of 
the SYM theory it effectively reduces to the Matrix quantum mechanics description of 
11-dimensional supergravity. Indeed, a new feature of matrix string theory (relative to 
standard light-cone string theory) is that via the electric flux of the gauge field, the string 
states can be adorned with an extra quantum number, identified with the D-particle charge 
. In a small g s expansion, these flux sectors energetically decouple, corresponding to the 
fact that D-particles can not be produced via perturbative string interactions. Nonetheless, 
electric flux can get created in the gauge theory: it is easy to see that electric flux creation 
is a simple one-loop effect that takes place whenever a virtual pair of charged particles gets 
created and annihilated, after forming a loop that winds one or more times around the a 
cylinder. 

In the following we will develop a new method for studying high energy scattering 
and D-pair production in Matrix string theory, which will be based on a semi-classical 
expansion from the SYM perspective. An important novelty of this method is that it 
applies to processes with arbitrary longitudinal momentum exchange. In the gauge theory 
language, this means that the transitions between the initial and final states that we will 
consider will involve a non-perturbative tunneling process in which an arbitrary number of 
eigenvalues get transferred between the two scattering states. Most previous calculations in 
Matrix theory relied on perturbative SYM corrections and thus were necessarily restricted 
to zero p + transfer^. 



2 In [16 1 Polchinski and Pouliot analyzed graviton scattering with non-zero M- momentum 
transfer in Matrix theory. In their case, the M-momentum was identified with the magnetic flux 
of the SYM gauge theory, and the corresponding instanton was a magnetic monopole. Here we 
will consider different kind of momentum transfer, namely of longitudinal momentum represented 
by the size N of the Matrix bound states, i.e. the number of D-particles in the original Matrix 
dictionary of @. This will require a different, less familiar type of instanton process. 
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Concretely, we will construct SYM saddle point configurations that will allow us to 
interpolate between ingoing matrix configurations of the form 

and outgoing configurations of the form 



xur) = ^ : J S ,* * r ) (i-3) 



+ 6)/ 3 
(fer-6)I 4 

where i$ are NixNi identity matrices, where all iVYs are different (but subject to the 
momentum constraint N1+N2 = N3+-N4). These mand out configurations each describe two 
widely separated gravitons with different light-cone momenta p^ = / R and transverse 
momenta pn\ , and with relative impact parameter b. 

The interpolating solutions that we will construct, essentially look like an appropriate 
matrix generalization of perturbative string worldsheets. The importance of these solu- 
tions is not entirely obvious, however, since a priori one would expect that the range of 
validity of the semi-classical Yang-Mills approximation has no overlap with that of pertur- 



bative string theory. Indeed, as emphasized in |T5[ the two regimes appear related via a 
strong/weak coupling duality. However, as we will argue in the following, even at small or 
moderate string coupling g s , at sufficiently high collision energies and/or impact parame- 
ters one enters a regime in which the semi-classical SYM methods may provide an accurate 
description of the scattering process. 

Just like string/M-theory, the 1+1 SYM model contains various length scales: (i) the 
circumference of the cylinder (set equal to 1), (ii) the scale set by the Yang-Mills coupling 

^YM 1/ 9yM 

Zym - 9s (1.4) 

(iii) the typical mass scale set by the Higgs expectation values of the SYM model. The 
latter length scale is inversely proportional to the impact parameter b of the string/M- 
theory scattering process: 

b 

Finally, (iv) there is also the length scale £e determined by the typical size of the SYM 
energy E, which is related to the relative space-time momenta via E ~ p 2 /N . 

The existence of these scales allows us to find small dimensionless ratios that may 
parameterize the strength of the SYM processes taking place at that scale. For example, 
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while g YM = l/g s defines the effective coupling of SYM processes that take place at the 
scale of the YM cylinder, we also have 



gfj(b) ~ 4/A 



1/6, 



[1.6) 



as the dimensionless coupling at the scale £5. Similarly, we can also associate an effective 
coupling g^J (E) with the scale set by the SYM energy E. This suggests the possibility 
that even if g s is small or of order 1, processes at these other 2D length scales can be 
accurately described by perturbative and/or semi-classical SYM methods. This will require 
however that we consider the limit of high collision energies and sufficiently large impact 
parameters.! 
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Fig. 2: The duality chain that illustrates the dual interpretations of the rank N 
and electric flux E in matrix string theory. 



The two types of processes that we will consider, high energy scattering with non- 
zero Ap_|_ and the D-pair production, may at first sight seem quite unrelated. However, 
there are several connections between these two types of processes. First of all, it is worth 
pointing out that in both cases the scattering process involves (depending on which duality 
frame one chooses) the transfer of D-particle charge and/or momentum between the two 
scattering particles. Indeed, the rank N started out as identified with D-particle charge, 

3 In this context it may be of relevance that in classical 10-dimensional DLCQ supergravity, 
the impact parameter b scales with the transverse relative momentum p via b ~ {^g 2 p 2 /N sinO^j 1 ^ 
with ./V the DLCQ p+-momentum. Hence, at least in this classical context, and for fixed scattering 
angles 6 and g s of order 1, the condition that b is large is automatically satisfied in limit of large 
p 2 > N. 
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and only after the duality it and the electric flux E are mapped onto each other under an 1 1- 
9 flip: i.e. the interchange of the 11-th and 9-th direction. (See fig 2.) Hence quantitative 
understanding both types of processes will have a direct bearing on the Lorentz invariance 
of the Matrix formalism. 

Another connection between the two calculations is related to a fundamental puzzle in 
the original calculation of 0, namely the apparently dominant contribution of arbitrarily 
high genus to scattering amplitudes. The saddle point trajectory at loop order G typically 
describes a process as depicted in fig. 3: two incoming strings, that are wound N = G + 1 
times, interact and then propagate as iV intermediate short strings. The iV strings then 
join together again, producing a final state of two different N times wound outgoing strings 
(see fig. 3). It was found in |I| that the contributions of these higher order interactions 
grows larger with the genus G. This instability appears to signal a fundamental breakdown 
of conventional string perturbation theory in the high energy regime. 




Fig. 3: This figure depicts a typical saddle point trajectory that contributes to 
the high energy scattering amplitude of fundamental strings, according to the per- 
turbative physical picture proposed in [jl| . 

On the other hand, the fragmented form of the intermediate state in fig. 3 gives a 
strong hint of some underlying non-perturbative structure that looks quite similar to that 
of the multi-D-particle bound state dynamics of Matrix theory. This suggests that the 
Matrix treatment may provide a rather natural stabilizing mechanism for a cutoff on the 
genus. Furthermore, our study will show that D-particle pair production becomes relevant 
at this cutoff - when the strings become maximally fragmented. This leads us to suspect 
a deeper relation between Gross and Mende's high energy, fixed angle scattering and the 
non-perturbative process of D-pair creation. 

We'll begin this paper with a quick review of some properties of Matrix strings. We 
then turn to fixed angle scattering. Section three will recall aspects of fixed angle scattering 
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in the traditional string framework, with particular emphasis on its description in light- 
cone gauge. This will be followed by a discussion of string interactions in the Matrix string 
formulation. In particular, we will find a local instanton solution in the two-dimensional 
theory that describes the splitting/joining interaction. Furthermore, the condition for this 
instanton to be matched to the incoming/outgoing states is precisely that we be working 
at the world-sheet moduli corresponding to the saddlepoint surfaces of 0. 

We then turn to D-pair production, which we consider both from the supergravity 
perspective as well as via a one loop Yang Mills calculation valid for arbitrary N. This 
is followed by a discussion of the ranges of validity of our calculations. In particular, we 
suggest that the ranges of validity of the two calculations overlap, and allow the picture 
we've mentioned in which they complement each other. This connection is further discussed 
in the final section, together with some other observations and speculations. An appendix 
contains discussion of the derivation of the Matrix string formalism following the approach 

of raj]. 



2. Matrix String Theory 

In this section we recall some basic features of the matrix string approach. This 
approach arises from DLCQ matrix theory quantized on a circle of "radius" R by com- 
pactifying on another circle, and interpreting this additional compactification as the route 
by which M-theory descends to type IIA string theory. For example, DO branes are states 
with non-zero momentum in this extra direction. 

2.1. Summary of Matrix Strings 

The Matrix string lagrangian is given in eq. ( |1 . 1|) . Its form can be derived from M- 
theory on S , by combining Sen and Seiberg's argument s [|17| , |T8|1 with the compactification 
prescription of fl9|] . (These arguments, and a more complete discussion of the relation 
among the parameters of the theory, are presented in the Appendix.) The basic dynamical 
variables of the theory are N x N hermitian matrices and include the 8 scalar fields X 
and the 8 fermion fields ^£ and i/j^. The Yang-Mills time variable r is related to target 
space time by r = tR in units where l s t = 1. 

The a direction corresponds to the T-dual of the compact direction (herein called x 9 ) 
of M-theory. The target-space coordinate in this circle direction is identified with the zero 
mode of the gauge field A a . The gauge equivalence under quantized shifts 

A a = A a + 27tm (2.1) 
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with m G Z translates via the identification of A a with the compact X 9 (see Appendix) 
into the periodic boundary condition along the M-theory circle. 

We can turn off the string interactions by considering only SYM configurations that 
describe widely separated strings, i.e. matrix configurations Xi such that all differences 
between their eigenvalues are large. In this situation all charged fields (relative to the 
Higgs scalar s) become very massive and effectively decouple from the dynamics. Hence in 
this limit all matrix fields take the form 



/ X[ \ 



X 1 



xi 



X N 



[2.2) 



and the corresponding effective SYM Hamiltonian reduces to the free field form 

N 

<J)da 

-9i 



1 - r r 1 i 
H =-22fda —El + (n£) 2 + {dvXlf + fermions (2.3) 
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where Ei denotes the (diagonal part of the) electric flux. 

For example, later in the paper we will consider matrices of the form 

pr + b ( I x /Ni 



where I\ and I2 are NixNi and N2 x N2 identity matrices. This corresponds (see 
Appendix) to two widely separated configurations with momenta fits = N^/R and 
P(i) = ±p/2, and with relative impact parameter b/2Ni + b/2N2- 

In the SYM/IIA string dictionary, the electric flux Ei gets translated into D-particle 
number. Indeed, in a canonical formalism the electric flux is conjugate to the (zero mode 
of the) gauge field A a , and thus represents the quantized momentum in the compactified 
direction. The quantum ground states of the SYM theory, corresponding to asymptotic 
particle states, are thus labeled by their transverse momentum p ± , their light-cone momenta 
p + and p~ , and their D-particle charge. 

Depending on their topology, the eigenvalue fields X- (a) in ( |2.2|) combine into one or 
more separate strings. For example, the trivial boundary condition 

Xi(a + l)=Xl(a) (2.5) 



8 



corresponds to a collection of elementary quanta (which may be thought of as minimal 
length strings), whereas the string of maximal length N is described by the periodicity 
condition 

X/(<r + l)=X/ +1 (<r), zG(l,...,iV). (2.6) 



We can write this condition as the N x N matrix equation 

X I (a + 1) = VX^^V' 1 



with V the cyclic permutation matrix on the iV eigenvalues, 



(2.7) 



V = 



V i 



i 





\ 



o / 



(2.1 



As a result of this cyclic boundary condition, we can glue the iV eigenvalue fields X-(o~) 
together into one single scalar field X 1 (a) defined on a long interval < a < N . Hence, 
when we expand this field in the usual way in oscillation modes, the frequency spacing 
between these modes is N times smaller than for a single valued field. 

In general the total matrix ( |2.2| ) will satisfy a periodicity condition of the form (|2.7| ) 



with V a block diagonal matrix consisting of (say) s blocks of order N/^ , such that each 



block can be taken of the form (2.8), and thus (as described above) defines a string of 
length Nuy As above, the space-time interpretation of this length is as the light-cone 
momentum 

N/.-\ 

(2.9) 



+ _ N H) 

P (i)- R 



This free string gas provides a good description of the SYM Hilbert space in the limit 
where all strings are far apart, i.e. if the eigenvalues in the matrix X are well separated 
between the different blocks. 

To avoid possible confusion later on, it is important to point out that the matrix V 
that specifies the periodicity condition on the eigenvalues is not equal to the Wilson line 
of the gauge field A a around the S 1 . This identification would arise only if we insist on 
minimizing the potential energy term (D a X) 2 with respect to A a . 

A related point is that for finite string coupling g s , and certainly in the large g s limit, 
the bound states with total light-cone momentum iV are clearly no longer necessarily 
described by means of a single long string. More generally, one would expect that the 
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bound state wave function will have support on more subtle bound state configurations 
consisting of several (up to N) short (er) strings. In other words, it seems reasonable 
to expect that for large g s the long strings will tend to "fractionate" into many smaller 
constituents. 

2.2. D-particles and Electric Flux Sectors 

As we have indicated, an important new feature of the matrix string formalism (rela- 
tive to standard light-cone string theory) is that via the electric flux, string states can also 
be adorned with a non-vanishing D-particle charge. In this subsection we will describe this 
correspondence in somewhat more detail. 

To add to this interpretation, let us first show that each separate string can carry only 
one type of electric flux. Consider again the single string with length N. Define the U(N) 
matrix U such that 

UV = VUe 2 ^ (2.10) 
with V as in (|2.8[) . Hence we can take 

/l \ 

2vri 



U 



2iri 

e n 



2(AT-l)7ri 

e n 



;2.ii) 



The SU(N) part of the electric flux in this sector is defined as 



U\ip e ) =exp(— -)|Ve) (2.12) 



N 

with e G Ztv and U the quantum operator that implements the constant gauge rotation 

(A, X) -> (UAU~\ UXU' 1 ). (2.13) 

Since diagonal matrices are inert under this gauge rotation, we conclude that the 
SU(N) part of the electric flux dynamically decouples from the diagonal configurations 
(|2.2|) that describe the separate freely propagating strings. Now recall that in U(N) SYM 
theory, the overall U(l) part of the electric flux is related to the SU(N) part e via 

tvE = e (mod N) . (2.14) 
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Supersymmetry ensures that the ground state in the SU(N) sector has zero energy even 
for e ^ 0. Hence the total ground state energy receives only a contribution from the overall 
U(l) flux. In the following we will thus identify e with the total U(l) electric flux. From 
the above description it is further clear that we can indeed turn on only one electric flux 
per long string, as is appropriate for its identification with D-particle charge. 



The energy (|2.3|) of the ground state in this electric flux sector is equal to 

ho = 4/ r < 2 - 15 > 

General ground state configurations 

|jV (i ),p«,eW) (2.16) 

of s separate strings of individual length ATW, transverse momenta and D-particle 
charge have a SYM energy equal to 



^o = E^iy^ ) ) 2 + (e W M) 2 l (2-17) 



i=l 



(recall this is defined with respect to the time r = tR) which, when rescaled by R, is the 
sum of the p_ light-cone momenta of the corresponding collection of string ground states 

£ -^\(P?) 2 + (M {i) ) 2 l (2.18) 
i=i 2 P+ 

In particular, we read off from ( |2.17| ) that the states with D-particle charge each have 



mass 



e (0 e (i) 



M« = — = — (2.19) 
9s R9 

in accordance with their identification as graviton states with non-zero KK momentum in 
the compact direction. 

2.3. IR limit and String Interactions 

In the IR limit, g s — > 0, the SYM dynamics effectively reduces to a free orbifold sigma 
model on the symmetric product space (R 8 ) N /Sn- The interacting theory for g s > arises 
by relaxing this IR limit. Correspondingly, one may view the interacting string theory as 
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obtained via a perturbation of the SAr-orbifold conformal field theory. In first order, this 
perturbation is described via a modification of the CFT Hamiltonian 



H = H + X J da V twist (2.20) 
Here V tw i s t is an appropriate CFT twist operator that generates simple transpositions of 



two string coordinates [[1511 . In terms of the bosonic twist-operators r and fermionic spin 



fields E it takes the form 



V twist = { tI ^i ® . (2.21) 



i<j 

This is a weight (f , f ) conformal field. The above twist field operator is an intertwiner 
between different topological sectors of the orbifold model that are related by a basic split- 
ting and joining interaction between two strings. Hence if we use the above Hamiltonian 
for computing scattering amplitudes via standard perturbation theory, we will indeed re- 



produce the conventional perturbation expansion of type IIA string theory [20]. This weak 



string coupling expansion is a strong coupling expansion from the SYM perspective. 



3. Fixed Angle Scattering of Strings 

High energy, fixed angle processes in superstring theory were first studied in detail 
from the point of view of conventional string perturbation theory by Gross and Mende . 
Central to their approach is the observation that in the limit of large external momenta, 
the Polyakov path integral at each given perturbative order is dominated by a finite num- 
ber of saddle point configurations. Furthermore, it was proposed that all these saddle 
points essentially describe the same preferred worldsheet trajectory, up to an overall factor 
depending on the loop order. 

In the subsequent sections we will find independent evidence from the point of view 
of matrix string theory that supports this physical picture. In this section we recall some 
of the main results of [1]]. In addition we will give a useful characterization of the Gross- 
Mende saddle points in terms of the light-cone gauge formulation of string perturbation 
theory. 
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Fig. 4: This figure indicates the kinematics of the transverse momenta pi. 



It will be useful to first recall a few geometric facts about the light-cone gauge formu- 
lation of string perturbation theoryE Consider a tree level string diagram that describes 
the scattering of four external massless particles with light-cone momenta pf = N^/R and 
transverse momenta pi. For definiteness, we will describe this process in the center of mass 
frame in the transversal direction 

Pi + P2 = 

(3.1) 

P3 + P4 = 

The four transversal momenta pi can all be chosen to lie within one given plane. We can 
thus write all pi as complex numbers. In addition, longitudinal momentum and energy 
conservation imply that 

N 1 +N 2 = N 3 + N 4 (3.2) 

|pi| 2 \P2\ 2 _ \P3\ 2 \P4\ 2 ( v 

N x + N 2 ~ N 3 + N 4 [6 - 6) 
For a given set of locations z% of the corresponding vertex operators, the classical 
location of the worldsheet is described by 

X+(z,z) = ij^iViloglz-^l 2 (3.4) 

i 

X(z,z) = - ^eiPi log \z - Zi\ 2 (3.5) 

i 

where ej = 1 for the incoming and —1 for the outgoing particles. Here (z,z) denotes a 
general conformal parameterization of the string world sheet. 



For a review of some salient features, see |21|. 
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In the light-cone gauge, one chooses a fixed world-sheet parameterization by identify- 
ing X + with the world-sheet time r 



X+(z,z)=t (3.6) 

which via (|3.4|) amounts to setting 

w = t + ia = — eiNi log(z - Zi) . (3.7) 

i 

The differential ui = dw is a specific globally defined holomorphic differential on the world- 
surface; existence and uniqueness of such a differential at arbitrary genus P^ , PT| generalizes 



the construction to higher loop amplitudes. Notice that (due to the branch cuts in the 
logarithm) the coordinate a in (|3.7|) is defined on an interval < a < (N\ + N 2 ). 

The light-cone coordinate system ( |3.7|) specifies a particular time-slicing of the string 
world-sheet. In this coordinate frame, there are therefore specific points on the world-sheet 
at which strings split or join. These interactions take place at zeros of u, that is critical 
points z = zq of the light-cone coordinate X + , at which 

dX+\ z=Zo =0. (3.8) 

Inserting the explicit form ( |3.4| ) for X + gives 

£-^V=0. (3.9) 

, Zq — Zi 
i=l 

In the specific case of the four-point scattering amplitude, this condition can be reduced 
to an equation relating the interaction point and the cross ratio 

A = (*l-*3)(*2-*4) (3 1Q) 

(21 - z 2 ){z 3 - z 4 ) 

via 

— + -[ = 4" (3 ' n) 

z Zq — l Zq A 

For given A, this is a quadratic equation for zq with in general two solutions Zq and Zq , 
representing the simple splitting and joining interaction respectively. 

Now we are ready to discuss the Gross-Mende saddle point. In the covariant for- 
mulation used in [0, it is characterized by the condition that it minimizes the Polyakov 
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action for the classical trajectory specified by (|3.4|) -( |3~5|) . This action takes the form of a 
two-dimensional "Coulomb energy" of four light-like "charges" given by the momenta pf. 

Ec = 7>^2 Pi-Pj log\zi- Zj\ 2 . (3.12) 

i<j 

Due to conformal invariance, this energy Ec depends on the locations Zi of the vertex- 
operators only by means of the cross ratio A. The variation of Ec with respect to A 
reads 

d x E c {\) = — ^— + ■ ( 3 - 13 ) 

The saddle point equation d\Ec = is solved by 

X = P -l^= t -. (3.14) 
Pi-Pn s 

This saddle point corresponds to a particular classical world-sheet trajectory which at high 
energies gives the dominant contribution to the scattering amplitude. 

For later reference, it will be useful to translate the above description of the GM saddle 
point into the light-cone gauge language. To begin with, in the complex parameterization 
for the pi, the Mandelstam parameters s and t are expressed as 



-2pi-P2 = (N 1 + N 2 



,2 



N^ 2 
\N 3 pi-N±p 3 \ 2 



(3.15) 

t = -2pi-p 3 



so that (|3.14|) takes the form 

_ N 2 \N3P1-Nml 

Together with ( p.ll| ), this saddle point specifies a particular set of locations for the two 
interaction points z^ of the light-cone string diagram. We now claim that this preferred 
location of the interaction points z = z^ is singled out by the requirement that, in the im- 
mediate neighborhood of z = z^ 1 , the transverse coordinate fields X(z) are (anti-)analytic 
functions of z 



d z X \ Z=Z0 =0 
djX \ Z=Z0 = . 
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(3.17) 



To verify this claim, let us compute the cross ratio A from (|3.17| ). The result should 
be equal to (|3.14| ). Inserting the solution (|3.5|) into (|3.17|) gives 

^ €iPi -=0 (3.18) 



In terms of the cross-ratio A defined in (|3.10| ) this reads 

Pi P3 



4(4-i) 4-a 



(3.19) 



where we used that P1+P2 = 0. When combined with the equation (|3.11| ), which relates A 
with the location of the interaction points zq, this equation can indeed be used to compute 
A in terms of the scattering data. If we subtract N 3 times (|3.19| ) from p 3 times (|3.11| ) , we 
obtain a linear equation for z , solved by 



+ _ N lP3 -N 3Pl 



Further, from ( |3.19|) we find that 



A = 4(l + ^(4-1)) . (3.21) 
\ Pi / 

After inserting ( |3.20| ) into (|3.21| ), it's a simple calculation to verify that the resulting 
expression for A indeed coincides with the high energy saddle point (|3.14| ). Note that for 
the saddle-point configuration, A is in fact real. The interaction points 4 an d are in 
this case each others complex conjugate. 

In [|l| , Gross and Mende propose the following attractive generalization of the saddle- 
point to higher orders in the string perturbation expansion. They assume that the domi- 
nant saddle points at genus G take the form of a G+l fold cover of the same four-punctured 
sphere as described above, branched over the four locations z% of the vertex operators. The 
resulting surfaces are known as Zn curves with N = G + 1. The classical trajectory of 
these higher order saddle points has the same shape as the tree level trajectory, but its size 
is N times smaller. (The intuitive reason is that they describe multiple wound strings, so 
that the effective string tension is iV times bigger than usual) . Correspondingly, since the 
different trajectories are weighted by the world sheet area, the higher order trajectories 
give contributions proportional to e~ Ec l N (with Ec given in ( |3.12| )). The higher genus 
contributions are thus quite strongly enhanced at high energy. We refer to jl] for a more 
detailed description. 
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It is worth pointing out that the structure of the curves and the corresponding 
space-time trajectories, as depicted in fig. 3, are quite reminiscent of the description of the 
"long string" boundary conditions in section 2. In our view, this (proposed) structure of 
the higher order interactions is one of several indications that the Gross-Mende approach 
to high energy string scattering may have a natural implementation in the Matrix string 
context. The above light-cone characterization of the GM saddle point in terms of the 
holomorphicity conditions (|3.17| ) will be critical in establishing this relation! 



4. Matrix String Interactions 

In this section we will prepare the ingredients for the semi-classical study of high 
energy scattering in the Matrix string framework. To begin with, we notice that the above 
light-cone gauge description of the dominant string world sheet trajectories can rather 
easily be put into a matrix form, as follows. Starting from equations ( |3.4]) and (|3.5|), 
we represent the classical string trajectory by means of a diagonal N x N matrix (with 
N = Ni + N2) by first writing the transversal coordinates X as a function of w defined in 
( P-7|) , and then "roll up" the spatial interval < a < N onto the short interval < a < 1. 
Concretely, we define the diagonal matrix elements of X(cr) via Xkk(&) = X(a + k), and 
in this way we indeed create matrix configurations that, away from the interaction times, 
satisfy the long string boundary condition ( |2.7|) and (p78|). 

These diagonal matrix configurations represent particular solutions to the SYM equa- 
tions of motion, that are regular everywhere except at the interaction points. If at some 
point in the (a, r) plane two eigenvalues Xj and Xj coincide, we enter a phase where 
locally the gauge symmetry is restored to U(2). In general we should thus expect that in 
this local region the semi-classical SYM solution will need to become truly non-abelian. 

It is readily seen that the diagonal matrix configurations constructed via the above 



procedure from the CFT solution (|3.5|) is not single valued around the interaction points. 
Instead, as explained in [IS], in going around the interaction point, the matrix X undergoes 
a simple transposition of the two degenerating eigenvalues. In the gauge theory language, 
this means that the diagonal CFT solution (|3.5|) in fact hides a delta-function Yang- 



Mills curvature at the interaction point, such that the infinitesimal Wilson line around it 



coincides with this permutation group element [23]. In this section we will describe how 



the Yang-Mills dynamics in fact smoothes out this singularity. 
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Fig. 5: The string interaction relating a one string to a two string state. This 
interaction occurs when two eigenvalues Xi and Xj coincide, we enter a phase 
where an unbroken U(2) symmetry is restored. 



Concretely, we will now exhibit a smooth and single-valued Yang-Mills configuration 
that describes the local splitting or joining of one or two matrix strings. Ultimately, we 
will be interested in obtaining global classical solutions to the SYM equations of motion 
that minimize the Yang-Mills action for given asymptotic conditions on the matrix fields 



4-1. SYM Solution near Interaction Point 

It seems reasonable to assume that, at least in the immediate neighborhood of the 
interaction point, these minimal action configurations of the SYM model are described by 
supersymmetric configurations. Hence, instead of trying to solve the full Yang-Mills equa- 
tions, we will restrict ourselves to the special class of solutions satisfying a dimensionally 
reduced version of the self-duality equations from four to two dimensions. We will choose 
to work with variables 



X, as written in eqns ( |1.2[ ) and ( |1.3j ) in the Introduction. 



X = -(X 1 + 1X 2 ) , X = -(X 1 - 1X 2 ) , 



(4.1) 



setting the remaining X l 's to zero. The self-duality conditions then become 



WW 








(4.2) 



. 



The above equations are most conveniently analyzed by writing 



A w (w,w) = —iGd w G 
Ayj(w,w) = i{dyjG 1 )G 



(4.3) 
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where G(w, w) denotes an element of the complexified (G 7^ G 1 ) gauge group. This 
parametrization of A a allows one to solve the second and third equation of (4.2), via 



X(w,w) = GX(w)G- 1 . (4.4) 



The first equation in ( |4.2|) then takes the following form 



with 



^(fid^fi- 1 ) = -^[nX(w)n-\ X(w)] (4.5) 
9 s 



O = GG. (4.6) 



Let us now look at the local neighborhood of an interaction point. For convenience, 
we choose coordinates such that it is located at w = 0. Since the interaction involves 
only two eigenvalues, it is sufficient to consider only the corresponding SU(2) part of the 
matrices. The matrix X, which parametrizes the local coordinate distance between the 
two interacting strings, can be chosen of the following form 



I(f)-±BVwr 3 (4.7) 
for some constant B. The ± indicates that the interaction point w = represents a square 



root branch point for the diagonal matrix X in (|4.7|) , which therefore is multi- valued. 

The diagonal matrix X(w), together with A = 0, represents a valid solution of the 
SYM equations ( f4.2| ) except at the interaction point, where analyticity fails. Therefore 
we'll look for a true solution of the form ( f4.4|) , where G — > 1 asymptotically far from 
w = 0. A helpful Ansatz for G{w,w) is 

G = e^ aTl (4.8) 

where for a(w,w) we choose a real function (so that G = G and Q = exp(ari)) that tends 
to zero far away from the interaction point. We now compute 

^>^^>-i 7-. /= cti -o-n 7-. /= ( cosh 2a — sinh2a \ , . n . 

QXfl 1 = B\fwe aTl Tze aTl = BVw . , n , n . (4.9) 

\ smh 2a — cosh 2a ' 



Hence 



2\B\ 2 \w\ sinh2a t x (4.10) 
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and thus we find that under the present Ansatz the equation of motion (O) reduces to 



d w c\ja = — o LB \w\ sinh2a 
9t 



(4.11) 



which is essentially the familiar sinh-Gordon equation. (It can be transformed to the exact 
sinh-Gordon equation after a ( mult i- valued) coordinate transformation w — > w = w 3 ^ 2 .) 

The boundary condition that we must impose on a(w,w) at w = follows from the 
requirement that the Yang-Mills configuration be regular. This condition is most easily 
understood in the gauge where X is single-valued near w = 0; in this gauge the YM 
curvature F„ 



should be a regular function at w = 0. The configuration ( |4.4| )-(p 
however, is (for single-valued and real a) multi-valued. We can make X single-valued by 
applying the singular gauge transformation 

X -> uxu~ l 

A w -> -iUD w U- x 

with gauge parameter 



(4.12) 



U = e ±l0T ^ 



with 6 



2 1 



log(w/w) the azimuthal angle around w = 0. In this gauge 



^ + 1 



Tl 



Using that d w l = n5^( 



w 



this gives 



— % 



(4.13) 



(4.14) 



n 



F - 

± WW 



(W 



(4.15) 



-iTi(d VJ d iS a± ^ (2) ( ? 

The regularity requirement at w = is therefore that d w d^ja ~ ^fj6^ 2 \w). We thus deduce 
that the solution to equation ( fl.ll|) that we want must satisfy the following asymptotic 
condition 

a(w,w) ~ =F- log \w\ + const. w^O (4.16) 

while at large distances from the interaction point a must tend to zero. 



Now let us write a = a(r) with r = 

ordinary non-linear differential equation 

1 



w 



. The equation of motion Q4.11|) reduces to the 



(d; + -d r )a = — |BMrsinh2a . 

r gj 

A numerical solution to this equation is depicted in fig. 6. For large r 
looks like 



a(r) 



| W |3/4 



exp 



8|B| 

3g s 



\w 



3/2 



(4.17) 
\w\ the solution 

(4.18) 
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Fig. 6: The numerical solution a(s) to (dg + ^d s )a = ssinh2a, as a function 
of the distance s from the interaction point. The initial condition for small s, 
that after integrating leads to the correct asymptotic behavior for large s, reads 
a (s) = -f logs + 0.0305070(1). 

5. High Energy Scattering of Matrix Strings 

The above matrix solution of the string interaction should be viewed as a local descrip- 
tion in the immediate neighborhood of the interaction point. In general, it must therefore 
be glued via an appropriate patching procedure into a complementary CFT type solution 
(e.g. as described in section 3) that matches with the asymptotic scattering data at the far 
past and future. The idea here is that (as we will see shortly) at sufficiently high collision 
energies, the size of the interaction regions are small compared to the rest of the matrix 
string world sheet. Hence, while the behavior ( f4.7|) provides the asymptotics for large \w\ 
at the UV scale of the matrix solution, it also provides the local boundary condition near 
the interaction point for the CFT solution for X that describes the IR part of the saddle 
point. 

The solution (|4.2|) that we described is not the most general SYM description of a 
string splitting and joining event, but rather the most symmetric one, with smallest action. 
This means therefore that there is a non-trivial matching condition on the corresponding 
matrix string world sheet: from (|4.7|) we see that we must require that the transverse 
string coordinates X behave holomorphically near the interaction point. Remarkably, this 
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is exactly the same condition as ( |3.17 ), which tells us that the shape of the string world 



sheet must be precisely that of the Gross- Mende saddle point! Therefore these solutions 
indeed seem appropriate to a YM generalization of the high-energy scattering of [JJ. In 
this section, we fill in a few more details of this connection. 

5.1. Evaluation of the classical action 

In order to estimate scattering amplitudes via the instanton processes, one must cal- 
culate the instanton action. The bosonic part of the SYM action (with only two X-fields 
non-vanishing) can be written as 

S = Jd 2 w j-s s 2 (i^ + A[x,X]) 2 + 4D w XD w Xj 



+ (j)(XDX + XDX - XDX - XDX) 



(5.1) 



and thus for the supersymmetric configuration that satisfy (|4.2| ) , the total classical action 
reduces to a boundary term 

S d = j)(XdX + XdX) (5.2) 

identical to the boundary term needed to glue the non-abelian matrix solution described in 
the previous section into the CFT type solution. Hence we claim that, in the limit that the 
matrix interaction points become sufficiently small, the SYM action for the above saddle 
point configurations coincides with the CFT action, i.e. for the case of a tree level string 
diagram it equals the "Coulomb energy" ( |3.12| ), where me must insert the saddle-point 
value for locations Zi. It is perhaps worth pointing out that this saddle point actions is 
indeed fully Lorentz invariant, as it should be. While this is not surprising once we've 
established the connection with the GM saddle point, it does seem to represent a rather 
non-trivial statement from the SYM point of view! 

More generally we see that from ( |5.1|) we can derive (as usual) an inequality, which 
suggests that whenever the interaction does not take place at a holomorphic point for 
the X-fields, the SYM action is always larger than the corresponding CFT action. This 
provides additional evidence for the conjecture that the above type of configurations indeed 
represent dominant saddle-points, that minimize the SYM action. 

Obviously, there exist a large number of CFT-type solutions for which X varies (anti-) 
holomorphically near all interaction points. In particular, there are the higher genus Zjq- 
curves of [1]]. In addition it is also possible to write down SYM solutions that describe 
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multiple string world sheets, but nonetheless still satisfy the appropriate boundary condi- 
tions, as specified in eqns ( |1.2|) and ( |1.3|) in the Introduction. Ideally, one would like to 
know which (sub-class) of these solutions provide the truly dominant contribution to the 
scattering amplitude. We will not attempt to answer this question here, and will restrict 
ourselves to some qualitative remarks in the concluding section. 

5.2. Minimal distance 

The parameter \B\ that governs the size of the interaction vertex, as seen in ( f4.18|) , 
can be straightforwardly determined in terms of the momenta of the external states. The 
coordinate system (w,w) on the Yang-Mills cylinder that we used in the analysis of the 
self-dual Yang-Mills equation ( |4.2|) , coincides with the light-cone coordinates defined in 
Q3.7Q . From this and Q4.7Q we immediately find 

,d|2 _ o \^ X (4) I 2 on 

A straightforward calculation then gives 

, B ,2 _ \PlP3-PlP3\(Nl+ N 2 ) 

It is interesting to note that for this solution, even though the eigenvalues of the com- 
plex coordinate matrix X indeed vanish at the interaction point, the full matrix coordinate 
X in fact does not! Instead, near w = it approaches the constant non-diagonalizable 
matrix 

X(w,w) ~ const. gl /3 B 2/s w^O (5.5) 

(The value of the overall constant can quite easily be determined numerically.) From this 
we read off that the minimal "distance" between the two interacting strings is in fact 
non-zero! Instead, we have 



drmn = y / tr(X(0)X(0)) ~ ^ 1/3 |S| 2/3 . (5.6) 

Although it is tempting to speculate (as indeed we will do in the concluding section), the 
precise physical significance of this result is as yet unclear to us. We do notice, however, 
that the typical world sheet size £i ns t of the matrix interaction region, as can be read off 
from ( f4.18| ), is naturally expressed in terms of this minimal relative distance as £i ns t = 

(Wl#l) 2/3 =9 8 /dmin- 
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5.3. Fluctuation Determinant 

In principle it should be possible to compute the one-loop determinant of the quan- 
tum fluctuations around these semi-classical configurations. An important motivation for 
performing such an analysis is to obtain a semi-classical estimate for the absolute strength 
of the splitting and joining interactions in matrix string theory. Duality symmetries of 
M-theory gives the precise prediction that this strength should be governed by the string 
coupling g s . To verify this, one would need to compare the SYM one-loop determinant 
with the Gross- Mende fluctuation determinant, coming from the gaussian integration over 
the Riemann surface moduli around the saddle-point. We leave this for future work. 

It seems even more worthwhile to look for true new physical effects that might arise 
from the one-loop corrections. Compared to the conventional perturbative string descrip- 
tion, the new degrees of freedom in matrix string theory are the charged components of 
the -Y-fields, as well as the extra gauge potential A a . These new degrees of freedom are 
non-perturbative from the string perspective, and their quantum fluctuations could thus 
potentially lead to new physics. As we will show in the next section, there is indeed such 
a new effect: the pair creation of D-particles. 



6. D-Particle Pair Production 

In this section we turn to the process of pair creation of D charge, which is in our 
description x 9 momentum or equivalently (under the Matrix string duality) electric flux. 
This can be viewed as a contribution to the fluctuations about the high-energy scattering 
processes of the preceding sections, or as a process worthy of interest in its own right in 
the context of graviton scattering. There are several viewpoints from which this can be 
investigated. In the limit where a; 9 decompactifies, this simply matches onto the standard 
supergravity calculation^?!]]. In fact, we can work backwards from this, using the method 



of images, and compute the amplitude at large finite Rg, in the special situation with 
source and probe particles, iVj ^> N 2 . We will discuss this calculation first. Alternately, 
one can study this process directly in the matrix string approach, and derive the pair- 
production rate via a one-loop Yang Mills calculation. This latter approach gives a leading 
order result valid for arbitrary N\ and N 2 , and also more readily makes connection with 
the other results of this paper. Furthermore, the Yang-Mills calculation also apparently 
extends beyond the region where supergravity is a valid approximation. 
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6.1. Supergravity calculation 

Consider 11-dimerisiorial supergravity compactified on S 1 xS 1 , where one S 1 a lightlike 
circle of "radius" R 

x~ =x~ + AnR (6.1) 

and the other S 1 denotes a space-like circle of radius Rg. As we've seen, in the M- 
theory /matrix string correspondence this second radius R9 is expressed as Rg = g s in 
string units. 

Consider in this set-up the scattering process of two massless particles of light-cone 
momenta pf = N^JR and transverse momenta pj. Let us first consider the probe situa- 
tion N2 <C N\. Then one can already get quite useful information about the scattering 
process from considering the classical gravitational force between the two particles. The 
boosted particle with p + = Ni/R produces via its stress-energy a non-trivial gravitational 
background, described by the generalized Aichelburg-Sexl shockwave geometry of the form 
23,^6 1 

ds 2 = —dx~ (dx + + f(r, a)dx~) + dx\ + g 2 s da 2 . (6.2) 

with 

Ar ' a> Y 2i? 2 (r 2 + g*(2nk + a) 2 ) 7 / 2 { } 

Here a denotes the coordinate distance from the gravitational source in the compact x 9 
direction, and the sum over k arises from the image points in this direction. 

The momentum four vector of the second massless particle moving in this background 
geometry will satisfy a dispersion relation of the form 

2p-(p+ + f(r, a)p~) =p 2 + e 2 /g 2 (6.4) 

where e denotes the quantized momentum in the xg direction. We can solve for the light- 
cone hamiltonian p~ of the particle and obtain 

P~ = 2^ <*/(!- Z ^Z^(P 2 + e 2 /9 2 ) ) - 1 } (6-5) 

Substituting p + = N2/R, (and rescaling the light-cone time by a factor of R) we can write 
this as 

H = H + H int (6.6) 



^{/MgW^)- 1 } 
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where 



Ho = ^- 2 ( P 2 + e 2 /g 2 s ) (6.7) 



and 

^"-^-4.7- — ^ + '-- (6 - 8) 

Hence the motion of the second particle in terms of the light-cone time x + looks like that of 
a particle with mass N2 moving in R 8 x S 1 under the influence of an interaction potential 
given by ( |6.8|) . 

From this description we can now quite easily extract a low energy prediction for the 
D-pair production rate. To this end, it is useful to rewrite the interaction Hamiltonian via 
a Poisson resummation as 

Hint ^ -^^g 2 (p 2 +e 2 /g 2 ) 2 ^M^a) I dTT 2 exp(-Tr 2 ) exp(-n 2 /4^ 2 T) (6.9) 

The n = 1 term in this series is the term that corresponds to changing the compact 
momentum by one unit, i.e. to D-charge production. Working to first order in perturbation 
theory, we can then compute the corresponding phase shift, using 

dTH tnt (b 2 +p 2 T 2 ) . (6.10) 

6.2. D-pair production via electric flux creation 

We now study this problem of D-charge creation in the Matrix string framework. More 
generally, we consider scattering states which asymptotically have momenta of the form 

= (p-,p,p 9 = n/R 9 ,p+ = N/R) . (6.11) 

These include both gravitons (n = 0) and DO-branes - or anti-branes - (n = ±1). The 
case of current interest begins with an initial state of two gravitons, and pair produces a 
D particle pair. This process is intrinsically non-perturbative from the point of view of 
string theory. It is also a process not accessible in the standard Matrix theory approach, 
where the anti-branes are boosted away to infinite energy. 

In principle (for example on a sufficiently large computer) it appears possible to calcu- 
late such amplitudes to arbitrary order in the coupling g = gyu = Vfi's, and calculate the 
D-pair production rate even for small g s . In this section we will work to leading non-trivial 
order (one-loop), and leave further calculations to other work. Similar calculations have 
been performed in the context of Matrix theory in E7 . 
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6.3. One-loop calculation, N = 2 

For simplicity we begin with the case where the incoming and outgoing particles all 
have N = 1. The next subsection will generalize to arbitrary N. The asymptotic states 
thus take the form fl2.4j) ; specifically, 



2 V -pr J ' 2 V° 

corresponding to two particles with center of mass momentum p and impact parameter b 
(measured in string units). It will also be useful work with a non-trivial gauge background 

+ °_ ; AW« M. (6.13) 



-a/2 - er/g 2 s J ' \0 -e 

The constant electric field corresponds to a non-zero D-charge for the incoming and out- 
going particles, with quantization 

e G Z . (6.14) 

The prototypical example of production of D-charge is in processes where this changes by 
one unit, 

AE=(l V (6.15) 



x _1 , 

Introducing the constant background potential a will help keep track of such changes. 

In the large string coupling/small Yang-Mills coupling limit, the leading contribu- 
tion to D-charge producing processes is easily computed via a one-loop super- Yang Mills 
calculation. Calculations at higher loop order then give subleading corrections in g = l/g s . 

Our starting point is the Yang-Mills action ( |1 . 1|) , although it will be simpler to begin 
with it written in its un-dimensionally reduced form in terms of the gauge field 

A M = (A»,gX i ), (6.16) 

with M = 0, ■ ■ • , 9, \i = 0, 9 and % = 1, ■ • •, 8. We will decompose the gauge field into 
background and fluctuation pieces, 

A M = A M + gA M . (6.17) 
The Feynman background gauge-fixed lagrangian is 

£ = -Tr ^Tr{Fl IN ) + ^Tr{D M A M f-i^pi^ (6.18) 
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where Dm = Qm + IAm- Using the decomposition (|S.17|) we find 
C = -Tr{ |^ + \{D M A N f + iF mn [Ami 

+ g^ + igD M A N [A M , A N ]- 9 -[A M , A n ] 



(6.19) 



The amplitude in question is given by 

A{a, e) = J VA^VX*ViPe lS (6.20) 

where the boundary conditions on the functional integral are chosen to correspond to the 
asymptotic behavior given in ( |6.12 ) , ( |6.13[ ) . 



If we write 



A M = ^A Ma <r a = ^A M+ <7 + + ^A M _a +^A M3 cr 3 , (6.21) 



with 

. a 1 ± ia 2 

then the couplings in ( |6.19|) include the standard charged minimal couplings of A + , A_, 
and to the U(l) field A^. The amplitude to create unit electric flux is therefore 
given by summing the loop contributions to ( |6.20| ) in which one of these charged particles 
circulates once around the a-direction; higher encirclings yield more flux. Therefore we 
need the contribution of the charged state windings to the one-loop amplitude. 

This immediately follows by reading off the spectrum from the second through fourth 
terms of ( B.19|) in the backgrounds ( |6.12| ) and ( |6.13|) . We begin by defining 

ft = g 2 p 2 + 4g 4 e 2 . (6.23) 
In the bosonic sector we find the neutral fields 



Jtf, i = !,-••, 8; m 2 = 



Al, m 2 = 

and the charged fields 

Xf , % = 2, • • • , 8 ; m 2 = r 2 = ftr 2 + a 2 + g 2 b 2 

\(gpA 9± -2eg 2 X l± ) ; m 2 = r 2 
V 

A 0± + -(2eg 2 A 9± +gpX 1± ) ; m 2 = r 2 + 2p 
V 

A 0± --(2eg' 2 A 9± +g P X 1± ) ; m 2 = r 2 - 2p. 
p 
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(6.24) 



(6.25) 



For the fermions, we have 16 massles uncharged states, 8 charged states with masses 
m 2 = r 2 + p, and 8 charged states with masses m 2 = r 2 — p, as in J27] . Finally, including 
the ghosts gives one complex, uncharged field C 3 with m 2 = and one complex charged 
field C ± with m 2 = r 2 . 

All of the charged fields are minimally coupled to the background field Ag = A a . At 
one loop level, we have 

Ai(a,e) = f n^ie iS(2)[$l] (6.26) 
J i 

where I labels the charged fields enumerated above (the uncharged contributions cancel), 
arbitrary winding is allowed, and where is the quadratic part of the action ( |6.19| ) 
including the coupling to through D. Working with phase shifts, we then have 



i5l = inAi = J2 ln H I V®ie iS(2H * l] 

I n J™ 

I n JO S Jn 



(6-27) 



Here we have used the functional integral representation in terms of the first-quantized 
trajectory cr M (s) = (r(s),cr(s)), n is the winding number about the cylinder, and Fj 
denotes fermion number of the field. 

For general winding n the functional integrals in ( |6.27|) are readily rewritten in terms 
of functional determinants. For example, with m 2 = r 2 we have 



I- 



VrVa e So ds ^ m2 /2-3 2 (p 2 t 2 +& 2 )/2-A^( s )] 

_ e -ina+in 2 /2S-ig 2 b 2 S/2^ p ^ ^ g\ 

where 



(6.28) 



A( P ,e >S )=de t -V 2 (^V (6 . 29) 



Combining such expressions and defining S = 2T then gives 



iS 1 (a^e) = e~ ina H ^^^Afe e, T) 

Jo 1 



6 — 2cos(2gpT) + 8cos(gpT) 

(6.30) 
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Recalling the quantization rule (|6.14|) , we see that as long as gp ^> 1 the electric background 
only contributes at higher order in g; neglecting this, the determinant is readily evaluated 
using 



deth(-d 2 s + g 2 p 2 ) =9P]J 



det 2 i{-d 2 ) 
The phase shift then becomes 



n=l 

2%T 



2mr 



2„2 



+ 9 P 



2sinh (gpT) 



(6.31) 



i6i(a,p,e) 



1 / i x - 
2V 



^ in 2 /4T-ig 2 b 2 T 



'T smh(gpT) 
-6 — 2cos(2gpT) + 8cos(gpT) 



(6.32) 



In the full amplitude A± = exp{z5i}, the coefficient of the term e ina is the amplitude 
to make a transition from a state with electric flux e to e + n: 



^i(a,e) = J2 e ~ ina ( e + n \ e ) 



(6.33) 



We have found that this is independent of e to order g 2 . The amplitude for a change 
by one unit of charge (e.g. two gravitons to DD pair), as well as the effective interaction 
hamiltonian, can be derived from these expressions in the range p b 2 . There the 
integrand in ( |6.32j ) can be expanded in pT to find 



iS 1 (a,p,e) « -^^Ij2 e ~ lUa J dTT 3 / 2 e™ 2 / 4T -^ b2T ; 



(6.34) 



the leading order DD production amplitude is just the coefficient of e~ ia in the series. 
From (|6.34| ) and fl6.10| ) we can also work backwards to extract the effective hamiltonian. 
We find 

1 



H 



int 



15 _3 4 

k 9 p 2^ 



8 Y l r2 + (a + 2nk)' 2 /g 2 ] 7/2 



(6.35) 



in agreement with (|6.8| ). 
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6.4- Generalization to arbitrary N 

The one-loop calculation of the preceding subsection is readily generalized to the case 
where the incoming and outgoing particles have arbitrary pn, or equivalently, N. In this 
case there are a variety of different boundary conditions that may be placed on the NxN 
blocks. Two that we will consider are the trivial boundary condition, 

X(a + 1) = X(a) , (6.36) 

and the single long string boundary condition, 

X(a + 1) = V~ 1 X(a)V (6.37) 

where V is given in (|2.8| ). 

In the case of two incoming states with momenta N%, we write 

X l = X l + X l (6.38) 

where X 1 is an (N\ + N2)x(Ni + N2) matrix. In particular, the background is taken to be 

^!_pr (I/Nx ) = pr 

2 I -I/N 2 J " 2 ' 



X 2 = -T D 
2 



D 

(6.39) 



where we have split the matrix into N\xN\ and N 2 x N2 blocks corresponding to two 
"clusters," and / represents the corresponding identity matrices. 

A useful decomposition of the fluctuations X % is in terms of the matrices 

r °' = ('o o)- r Mo £)■ < 6 - 4o > 

where t ai are hermitian generators of SU(Ni);T" 1012 ,T^ ia2 , which have matrix elements 

(6.41) 

(T 2 1 a ) /3l/ g 2 - V / 25 A r 1+Ql/ 3 1 5 a2/ 3 2 



32 • 



and Tp: 



-X-D rp Xqi rpg x 1 rpg 2 . ^+ai02ma 1 a 2 , _ " I " . .- 



_ 2 D 2 2 2 + 2"' y • > 
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Following the preceding subsection (and working with e = for simplicity) we find 
that the charged states now have an extra N1N2 in their multiplicities, and have masses 





X i 


m 2 = 




A 9 

±OL\OL2 ' 


m 2 = 


A , 


4. x 1 


m 2 = 


A , 


- X 1 


m 2 = 



Av 2 

f 

r 2 gp 
Av 2 v 



(6.43) 



r 



Au z 



gp 



V 



where 



111 

v ~ + N2 



(6.44) 



Likewise, the charged fermions and ghosts have masses as in the N = 2 case with the 
trivial rescalings to 

(6.45) 



p r b 
P = — i = — 



Therefore, in the case of trivial boundary conditions the amplitude (and hamiltonian) 
is exactly as computed in (|6.32| ) with the only difference being multiplication by N1N2 and 
replacement of p and b as in (|6.45|) . Note that 2p = + ]y- is simply relative velocity 
of the two clusters, and b = + jfe) is precisely the impact parameter between the 

clusters. 

In the case of long-string boundary conditions, this result is modified. Now 



(6.46) 



and in particular the charged off-diagonal blocks satisfy twisted boundary conditions 

(6.47) 

The matrices V can be diagonalized by working on basis vectors 



X+(a + l) =V l - 1 X+(a)V 2 
X-(a + l) = V 2 ~ 1 X-(a)V 1 



1 



1 

\ k 
\ 2k 



A = e 2ni / N keZ 



(6.48) 
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and in this basis simply give phases X k . Thus the amplitude (|6.32|) is modified to 



TV, 



iS 1 = - J— ^ e~ in[27T(cei/Nl ~ a2/N2)+a] [ ^ c in-/-iT->!flrT 



2 V 2yr c^i n J o T VTsm(gpT) 

—6 — 2 cos(2gpT) + 8 cos(gpT) 



(6.49) 



and the interaction hamiltonian takes the form 



Hint W ^ ^ e - l n[2.(a 1 /TV 1 -a 2 /Ar 2 )+a] ^ /4T-z, 2 T (g 5Q) 



n,ai,a 2 



For non-zero n, the supergravity correspondence no longer holds when N > 2: the 
matrix string then yields a different result. In ( |6.50| ) the expression in the summation only 
gives a non-zero contribution when the integer n is a multiple of both Ni and N2. Hence 
we see that the minimal exchanged D-particle number between two long strings of length 
Ni and N2 must be proportional to N1N2 (if the lengths are relatively prime), else the 
amplitude will be simply zero. 

This leads us to the conclusion that the long strings do not give an effective means 
of creating D-particles. For two strings to create a minimally charge D-pair, the strings 
apparently first need to each emit a minimal length string such that two short strings of 
both collections can exchange a single D-particle. In the SYM language, this last process 
is effectively an SU(2) process, where correspondence with 11-dimensional supergravity 
is found. It is important to note, however, that the electric flux thus created must then 
subsequently spread out over the complete U(Ni) gauge group, since otherwise it would 
not carry the SYM energy appropriate for the massive D-particle with p + = Ni/R. 

Furthermore, in the sector with a fixed p + momentum, we now have an improved 
idea of what state contributes most to D-charge production: it is the state with trivial 
boundary conditions, ( |6.36| ) , corresponding to a collection of minimal length strings. Since 
this is the state that yields amplitudes agreeing with low-energy supergravity in the limit 
g — > 0, it is apparently this state (or a bound version of it when finite g effects are taken 
into account) that dominates the wavefunction of the graviton in the small g region, rather 
than the state with the long string boundary conditions (|6.37| ). 
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7. Ranges of validity 

In this section we will give a preliminary discussion of the relevant scales and ranges of 
validity of the calculations of the preceding sections. This analysis is preliminary in that 
the systematics of the perturbation theory for the Yang-Mills lagrangian ( |1 . 1| ) has not 
been performed at the level of that for pure Matrix theory [p6| and additional subtleties 
are possible. We leave such analysis for future work. For simplicity we will consider the 
case where the p + momentum of the two incoming states are comparable, N± ~ N2 ~ N. 
Our arguments readily generalize to the probe situation N\ ^> N2 

We begin by considering the expansion of the action about a classical background as 
in ( |6.19|) ; such a treatment is relevant both for corrections to the saddlepoint solutions of 
section four as well as for the systematic treatment of D-pair production. 

This expansion is governed by the Yang-Mills coupling gyM, and naively one expects 
the condition gyM *C 1 for corrections to be small. However, as mentioned in the intro- 
duction, the Yang-Mills coupling is scale dependent and one expects the relevant scale to 
be set by the physics one is considering. 

For example, in the scattering with background ( |6.12|) , loops of the charged, massive 
states of the YM theory play a central role. One either has a loop localized on the cylinder, 
whose calculation leads to the 0{v 4 /r 6 ) supergravity potential, or the loop can encircle the 
cylinder leading to the D-pair production that we have computed. These massive states 
receive masses of minimum size b/g s through the Higgs effect, setting the length scale 
it, ~ g s /b. At this scale, we expect the relevant dimensionless parameter to be 

g YM lb ^ ^ • (7.1) 
Smallness of this parameter thus requires 

6»1. (7.2) 

For the case of pair creation, there is another requirement arising from the condition 
that the backreaction due to the created electric field be small. One way of stating this is 
to require that the YM energy be large as compared to the energy stored in the electric 
field, 

P 2 » 9ym ■ (7-3) 
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Finally, in the case of the string interactions of section four, we see from ( f4.18|) that 
the relevant scale is set by the parameter \B\, and is given by 

n 2/3 / 2sr \ V 3 

9s V I 9s N x 



At this scale the dimensionless coupling is given by 



N 



9YM^inst - o ■ a • (7.5) 

Another condition to apply the methods of section four is that the size of the instanton be 
small as compared to the size of the cylinder, li ns t <C 1, or 

P 2 » ^ • (7-6) 
suit/ 

It is certainly possible to simultaneously satisfy the conditions (|7.2|) , (|7.3|), and (|7.6|), 
as well as the more stringent condition gyM ^ 1; f° r finite N and large s ~ p 2 . If all 
important corrections are governed by expansions in the parameters of and (|7.5|) , 
then it appears possible to even push the calculations into the range g s <l. 

A more complete analysis can be performed in the large g s (large Rg) case in the 
restricted energy range 

— < E < g s . (7.7) 

9s 

The lower bound corresponds to the energy threshold to create D charge, and the upper 
bound is the energy to create winding states wrapping x 9 . In between these bounds the 
theory can be effectively described by Matrix theory DLCQ quantized in 10 dimensions. 
As explained in |26[], the Matrix expansion is an expansion in terms of the form 



N\ L ( v 2 



(7.8) 



where L counts loops. The terms with L = n are readily identified with terms in the 
corresponding supergravity expansion, and the small parameter justifying this expansion 

Nv 2 , s 

^ « 1 . (7.9) 

This has a simple physical interpretation, which is easily seen by estimating the net trans- 
verse momentum transfer due to the potential 
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this gives 

N 2 v 3 

Apx = w . (7.11) 

The condition ( [Tip is then easily seen to be Ap± <C p, or equivalently 9 <C 1 where 6> is 
the scattering angle. Combining this with the threshold condition ( |7.7| ) then yields 

<? s p » 1 , (7.12) 

in correspondence with (777). 



Expansion terms with n > L are then suppressed for 

p 2 < iV 2 r 4 , (7.13) 

and terms with n < L for 

r > (JVsO 1 / 3 . (7.14) 

It is unclear whether the latter condition is strictly necessary; the first term in this expan- 
sion vanishes pT| , |2"6[j , and the other terms have been conjectured to vanish in ||28|| . 



To better understand these conditions, we convert them into statements relating the 
Mandelstam parameter s ~ p 2 and the angle 9. It is easily seen that condition ( |7.13| ) 
becomes 

* « N 2 (j) (7.15) 

and the condition ( |7.14j ) becomes 

s » g 7 J 3 N 10 / 3 9 . (7.16) 

Comparing (|7.6|) with (|7.15|) and (|7.16|) , we see that within the energy range given by 
( [7.7|) the instanton and D-particle production calculations are not obviously simultaneously 
valid. However, outside of this range, we appeal to the preceding (less rigorous) analysis 
which suggests that these calculations are indeed simultaneously valid at large s, and may 
even be extendable to g s <l. It is partly on this basis that we will, in the next section, 
consider the consequences of combining these two calculations. 



36 



8. Discussion and Conclusions 



We begin this section by recalling several observations from our preceding discussion. 
The first is that, as pointed out in section 6.4, string scattering only efficiently produces 
D charge if the strings break off at least one minimal length string. Furthermore, sections 
four and five discussed saddlepoint configurations that are expected to make important 
contributions to high energy, fixed angle scattering. Combining these yields a picture of 
how the important non-perturbative process of D-charge production can arise in high- 
energy string scattering. 

The analysis of Gross and Mende|lj] found saddlepoints believed to dominate scattering 
at high energy. These saddlepoints have a common structure at arbitrary genus, and the 
contributions of these saddlepoints grows with the genus suggesting the relevance of non- 
perturbative effects. We have found a new version of their analysis in which a mechanism 
appears that can cut off this growth. The cutoff originates from the minimal string length, 
which is in our language the minimal p + . String fragmentation is stopped when the string 
breaks into the maximal number of minimal-length stringsJl 

It is precisely in the context of minimal length string scattering that we have found 
that D-charge pair production can become an important effect. We therefore have a very 
nice picture in which the instantons of section four and five lead to maximal fragmentation 
of the strings, and this is followed by the production of D-charge via the process of section 
six. Here we expect that the size ( f4.18| ) of the instanton, as well as the corresponding 



minimal distance ( p.6[ ), may be an important ingredient in determining the size of both 
these effects. 

From the stringy viewpoint this is an intrinsically non-perturbative process. This is 
suggestive that there is in fact a basic connection between these two processes, and in 
particular that the non-perturbative production of D-charge is an important correction to 
the high-energy scattering analysis of pp. While we believe that, by combining the various 
ingredients presented in this paper, it may be possible to obtain definite quantitative 
estimates of these corrections, we leave further analysis of this connection for future work. 

Next we turn to several other observations and connections. 



5 In this sense, the matrix string formalism behaves very similar to the discretized models of 
string theory, advocated in particular by C. Thorn. We thank S. Shenker for drawing our attention 
to this similarity. 
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First, recall that Banks and Susskind [^9| previously considered the DD system in 
the context of perturbative string theory. There they found an instability with unknown 
outcome. In the present framework we have been able to treat the same system analytically, 
at least in the large g s limit, without signs of pathology. In principle, the Matrix string 
calculations appear to extend to arbitrary g s . One might hope that some extrapolation of 
our approach could shed further light on the discussion of [p9| . 

It is an interesting conceptual question under which circumstances one needs to include 
the virtual effects of D-particles propagating in loops. Although in the literal sense of an 
expansion about g s = they do not contribute, since they have infinite mass there, there 
is clearly a strong sense in which D-particles can be found in intermediate states when g s is 
finite. Indeed, intermediate states with D-charge are distinguished from other intermediate 
states only by the presence of electric flux, and there is no apparent reason why these should 
be suppressed at finite g s . In fact, looking at the results of section 6, leads one to suspect 
that it may be possible to extend the matrix string interactions as summarized in section 
2.3 to include the possibility of electric flux "pair" creation. The eleven dimensional 
symmetry of M-theory, in particular, suggests as a possible generalization of the DVV 
string-interaction vertex, an expression of the form Vi n t = Vt w i s tS(Ai2) (with A12 the 
difference between the U(l) gauge fields on the two strings that are created). With this 
choice of vertex, the couplings between string and all n-D-particle bound states are all of 
the same strength. This would suggest that there may possibly exist a systematic semi- 
classical expansion in string theory - generalizing the standard perturbative expansion - in 
which the D-particle-loop contributions play the role of instanton-like corrections. Indeed, 
in other recent studies of non-perturbative contributions to string scattering amplitudes 
30| - |33[1 it was suggested that D-particle loops are related via T-duality to D-instanton 



contributions in IIB string theory. It clearly would be interesting to see if the suggestive 
formulas obtained in these works can possibly be reproduced via the Matrix string methods 
developed in this paper. 

To conclude, we have succeeded in using the Matrix String approach to begin an 
investigation of aspects of high energy string scattering, and in particular to begin to 
explore the role of important non-perturbative (from the string viewpoint) processes such 
as D-charge production. Further investigation along these lines is expected to unravel 
a rich structure at substringy scales, and may shed further light on the short distance 
structure and fundamental degrees of freedom and dynamics of M-theory. 
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Appendix A. The Matrix String Approach at Finite N. 

Most of our calculations in the main text are performed for type IIA string theory in 
the DLCQ approach, in other words for the IIA string compactified on a lightlike circle of 
"radius" R 

x~ = x~ +4ttR (A.l) 

This means that from the M-theory point of view we have compactified on S l x S , 
with the second compactification on a circle whose radius Rq is expressed in terms of the 
corresponding string length as Rg = g s £ s . Through most of the paper, we work in string 
units l s = 1, so that Rg = g s . 

In this Appendix we will rederive the 1+1 dimensional U(N) SYM lagrangian as the 
appropriate Matrix formulation of M-theory in this discrete light-cone limit. This is useful 
in understanding the precise scalings and relations between parameters. We will follow 
closely the reasoning of |T8| , pr7| , p!9| . 



To begin with, consider M-theory compactified on a circle with radius Rn <C l p i, 
where l v i is the eleven-dimensional Planck scale. Suppose we restrict attention to the 
sector with momentum pn = N/Rn. The lowest excitations in this sector occur at energies 
E ~ RuM^, and are described by the action (here for simplicity we omit fermions) 0,[5|,[| 

S = I 2^ (^ + \ Rlx ^ 1 [X ' X]2 ) • (A - 2) 
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X is an N x N matrix; the diagonal terms in X correspond to the coordinates of the branes, 
with kinetic energies MdqX 2 /2 = X 2 /2Rn and the off-diagonal components describe 
the creation of straight strings stretched between the branes, with characteristic energies 
~ \xi — Xj\RnMp L . This action neglects string oscillations and higher energy excitations 
(Planck scale excitations, brane creation). The lowest of these begin at the string scale 
(defined with respect to the compactification on Rn) 



M s = y/RuM* . (A.3) 
Next we compactify on a circle of radius J?g, through the identification Xg ~ Xg + 



2tvRq. In the matrix context, Taylor |nj has argued that this is most easily described by 
passing to the dual circle of radius R' 9 = \j2-KRgM 2 . The Xg variables dualize to a gauge 
field through the identification 

A x =M 2 J2 einX/R ' 9X 0n (A.4) 

n 

and the wrapped strings combine into a new definition of the remaining fields, 

X* = J2e inx/R ' a X % 0n , (A.5) 

n 

with the resulting action 

dt r 27TR '° dx „ f~2 1 so / „ . r . _x2 i 



rdt r* 9 J^ Tr ( x \ ^La 2 x - (a x x + i[A x ,x}) 2 + l-Mfix ,x} 2 \ 

J 2R U J 2nR' 9 \ M s 4 x V / 2 J J 



(A.6) 



It's then convenient to redefine the spatial coordinate, 



x = 2nR' 9 a ; A x = , (A.7) 



2tvR' 9 



giving 



S = J ^^^aTr|i%^-^(I) CT X) 2 + i^[X,X] 2 | . (A.8) 

The action ( |A.8| ) so far corresponds to the low-energy action of the system of N DO 
branes on the circle of radius -Rn- The next step is to take the limit 

R u -»• 0,M p z -> oo (A.9) 
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holding 



M^Rn = MiR (A.10) 



fixed, following |H|,|T7|]. This gives the infinite momentum limit, p\\ — > oo. At the same 



time, one defines a rescaled X and Rg, 

X = ^X ; R 9 = ^Rg . (A.ll) 

M pl M pl 1 1 

In the limit (|A.9| ), ([A.10| ) all the higher excitations mentioned above (string, Planck, and 
DO) decouple and the lagrangian ( |A.8| ) becomes exact. Written in terms of the rescaled 
variables, 

J JO I pi pi 



Eq. ( |A. 12 ) therefore describes DLCQ IIA string theory lightlike compactified on a 



circle of "radius" R. The physical string coupling and string scale are defined by the usual 
relations 

#9 = 9sht , Ipi = gl /S lst ■ (A.13) 
Finally, as a matter of convenience it is easier to work in terms of a rescaled time 

r = tR/l 2 st , (A.14) 

and to redefine 

X -> l st X , (A.15) 

corresponding to measuring X in string units. With these redefinitions the action indeed 
takes the form as given in the introduction. In the main text we will take l s t = 1. 

In order to describe physical scattering processes, we will also need to parameterize 
the external states. We will separate off the two compact momenta from the transverse 
momentum p, and write 

p» = (p-,ftpg = n/R 9 ,p+ = N/R) . (A. 16) 

For example, when we compute scattering of two gravitons into two DO particles (now 
defined with respect to the Rg compactification) , initially n = and finally n = ±1. 
Notice that scaling the momentum according to ( |A.11|) , finite p corresponds to infinite 



41 



p = M p ip/M p i. However, the corresponding energy scale in the original variables of (|A.8| ) 
is E ~ p 2 Rn = R 2 p 2 which remains finite as the other excitations decouple in the limit 
Rn — ► 0. In the case of two particles in the center of mass frame each with one unit of 
DLCQ momentum, and with relative momentum p. the asymptotic state is described by 

X = l( m n +i ° A- (A.17) 



2 V -Rpt - b 
In the rescaled units defined by (|A.14|) , ( [A.15Q , and l si = 1, this becomes 

x=\(fr+ % y A , (A.i8) 



2 V -pr - b 

with the understanding that the external momentum and impact parameter are measured 
in string units. 
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